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Abstract 


Currently, ionomer membranes are used in a variety of specialized applications. Such applications include, but are not limited to, dialysis, 
electrolysis, membrane separators, reaction catalysts and the most promising application: polymer electrolyte fuel cells. Although their use is 
widespread, significant gaps in understanding the mechanical behavior of these materials still remain. Many ionomer membranes change their 
structure, and in turn, their mechanical properties in response to applied thermal and moisture conditions that are functions of position. It has been 
observed that constrained materials subjected to changing environmental conditions can exhibit unusual behavior, e.g., in some cases, mechanical 
failure is seen in the absence of external applied mechanical loads. This condition is especially important in polymer membranes (specifically 
Nafion®) used in polymer electrolyte membrane (PEM) fuel cells and is the major motivation of the present work. Laboratory characterization 
has been conducted to determine the mechanical properties of a proton exchange membrane with respect to temperature and relative humidity. 
Data recovered in these tests along with properties from literature have been used in finite element models to predict the behavior of membranes 
used in certain applications and geometries. The overall goal of this investigation was to characterize the mechanical response of ionomer 
membranes in in-plane constraint configurations subjected to variable hygro-thermal environments. Expansion/contraction mechanical response of 
the constrained membrane as a result of change in hydration and temperature is studied in uniform and non-uniform geometries and environments. 
With this information, mechanical failure modes can be analyzed which is necessary for durability modeling and life prediction. The present work 
concentrates on defining and understanding the basic mechanical behavior of ionomeric membranes clamped in a rigid frame, and subjected to 
changes in temperature and humidification. 
© 2007 Elsevier B.V. All rights reserved. 
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1. Introduction 


Polymer electrolyte membrane (PEM) fuel cells, which uti- 
lize hydrogen to generate electricity through electrochemical 
reactions, may be an environmental-friendly replacement for 
current internal combustion engines in transportation (and other) 
applications. Durability and cost are the major barriers to the 
commercialization of PEM fuel cells. The DOE durability and 
cost targets (2015) set for an 80 kW (net) transportation fuel cell 
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stack (excluding balance of plant BOP) are a lifetime of 5000 h 
over a wide range of temperature and a per-kilowatt cost of 15 
dollars or less (which corresponds to 20 dollars per square meter 
for polymer electrolyte membranes) [1]. Membrane failure is 
one of the major stack life-limiting factors, which reduces the 
durability to the current 2000 h status (2006). To meet the DOE 
durability targets, fundamental understanding of the degradation 
and mechanical failure of ionomer membranes in the PEMFC 
stacks is necessary to help mitigate this type of failure. In order 
to establish a mechanistic durability and life prediction model, 
three steps have to be taken: (1) investigate the mechanical 
behavior of ionomer membrane; (2) establish accurate numerical 
stress—strain constitutive modeling methods for PEMFC stacks; 
(3) identify material’s degradation and failure mechanisms and 
kinetics for use in a mechanistic-based life prediction model. 
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This work will focus on the second part in conjunction with the 
previous investigations performed by the authors to tackle the 
outlined steps. 

The most commonly used membrane in PEM fuel cells is 
a perfluorosulfonic acid (PFSA) membrane. Among the PFSA 
polymers used in fuel cells the most favored one has been the 
sulfonated tetrafluoroethylene copolymer with the trade name 
Nafion®, developed and manufactured by DuPont™ in the 
1960s. It is the first of a class of synthetic polymers with 
ionic properties, which are called ionomers [2]. Nafion® has 
a very complex structure which has intrigued scientists for 
decades especially in the recent years when new applications 
in fuel cells appeared to be a solution to energy problems. 
Although there are a variety of models to describe Nafion®’s 
structure, it is still not very clear how it changes in different 
environments. Nafion® consists of a hydrophobic polytetraflu- 
oroethylene (PTFE)-like backbone and pendent chains with 
sulfonated (SO37) end groups. Under humidified conditions, 
the hydrophilic end groups segregate into nano-sized clusters, 
which imbibe water and cause the swelling of the ionomer. The 
phase-separated morphology and the ionic conducting mech- 
anisms of Nafion® have been extensively studied and more 
research is underway to understand the micro/nano behavior 
of perfluorosulphonated ionomers [3-15]. These water clusters 
form interconnected acidic domains, which are also thought to 
be the major path for proton transport. Both the swelling and the 
proton conductivity increase with the level of hydration [16-20]. 
The mechanical behavior of Nafion® is also known to be a strong 
function of hydration [21-26]. 

Membrane failure is believed to be the result of combined 
chemical and mechanical effects acting together [27,28]. Varia- 
tions in temperature and humidity during operation cause cyclic 
stresses and strains (mechanical loading) in the membrane as 
well as the MEA and is considered to be the mechanical fail- 
ure driving force in fuel cell applications [21,29-31]. Reactant 
gas cross over, hydrogen peroxide formation and movement, 
and cationic contaminants are all believed to be major factors 
contributing to the chemical decomposition of polymer elec- 
trolyte membranes. Huang et al. have observed localized attack 
of Nafion® under open circuit voltage (OCV) degradation lead- 
ing to crossover failure. This localized membrane attack begins 
on the cathode-side and is rapid in dry and hot conditions [28]. 
When exacerbated by mechanical effects, e.g., relative humid- 
ity cycling, it leads quickly to crossover failure as numerous 
pores/cracks propagate across the MEA and pinholes will be 
formed [20,27,32,33]. These pinholes cause gas crossover and 
localized heating of the membrane and can eventually lead to 
catastrophic failure of the fuel cell stacks. 

While chemical degradation of prefluorinated membranes has 
been investigated and reported extensively in literature [34—43], 
there has been little work published on mechanical degradation 
of the membrane. Cyclic hydration of the membrane/MEA dur- 
ing the operation cycles (start/shut down) of the fuel cell and 
hence mechanical fatigue of the membrane electrode assem- 
bly has been shown to cause mechanical degradation of the 
membrane [29]. Investigating the mechanical response of the 
membrane subjected to fuel cell cycles requires studying and 


modeling of the stress—strain behavior of membranes and MEAs. 
Weber and Newman were the first to study the stress associated 
with constraint in a simple one-dimensional model [44]. A finite 
element model has been used by Tang et al. to investigate the 
hygro-thermal stresses in a cell assembly caused by tempera- 
ture and relative humidity change with an assumption of linear 
elastic behavior for the membrane [31,45]. Later, Kusoglu et 
al. modified the work to simulate the stress evolution during a 
simplified operating cycle using a linear elastic-perfectly plas- 
tic material model. They have shown that compressive, plastic 
deformation occurs during the hygro-thermal loading, resulting 
in tensile residual stresses after unloading [30]. 

It has been shown that stress relaxation of Nafion® is highly 
influenced by introducing water in the structure; thus, any visco- 
elastic/plastic model for Nafion® should consider hydration 
effects in the stress model [46]. Benziger and Chia proposed that 
membrane swelling and relaxation processes work as an interfa- 
cial contact switch between the membrane and the catalyst layer. 
In this case stress relaxation causes the extruded membrane into 
the electrode to slowly relax and stretch out in the plane of mem- 
brane [47]. The importance of the relaxation and also the existing 
time delays during hydration/dehydration of the membrane in 
fuel cells reveals the significance of adding viscous contribution 
to the stress model. A visco-elastic stress model has been devel- 
oped by Lai et al. and stress predictions upon implementing 
the model using the data for Nafion® NR111 have been vali- 
dated with stress measurements from a relaxation test performed 
at small initial strain (3%) in the linear elastic region. They 
have also treated Nafion® as a rheologically simple material and 
applied time—temperature—moisture superposition principles to 
generate master curves for creep compliance of Nafion® [29]. 
In the authors’ previous studies, the properties obtained for a 
wide range of temperature and hydrations have been applied in 
a finite element model to study through-thickness stress—strain 
behavior of an MEA as a result of hydration and temperature 
cycling during MEA fabrication and/or unitization with GDLs. 
Geometric irregularities have been introduced to induce high 
local stress—strain concentrations, and as potential sources for 
crazing and crack-like defects leading to final mechanical failure 
of the MEA [21]. The same elastic—plastic material model has 
also been used for in-plane numerical modeling of RH-induced 
strain in an MEA in a constraint configuration [28]. In the present 
paper, a two-layer visco-plastic model is proposed to predict 
strain rate dependence and relaxation behavior of a membrane, 
and will also be discussed further in a future work. 

This paper focuses on a mechanical modeling scheme of 
the membrane in order to simulate the stress—strain response 
of the membrane in the constrained frame. This in-plane con- 
straint configuration roughly mimics the clamped membrane in 
the fuel cell. Changes in hydration and temperature will be the 
driving forces, and hydration will be shown to have a bigger 
effect. It is important to understand that without any exter- 
nal applied force, the membrane can be weakened primarily 
by chemical reactions and by cyclic hydration, and can subse- 
quently fail mechanically in the constrained geometry. Thus, 
hydration (and to lesser extent temperature) acts as the driving 
force for mechanical breach and formation of pinhole/cracks 


368 R. Solasi et al. / Journal of Power Sources 167 (2007) 366-377 


in the framed membrane. This article is outlined as follows. 
Mechanical properties of Nafion® are presented as a function of 
relative humidity and temperature. Then those material curves 
have been incorporated into the commercial finite element code 
ABAQUS (HKS). Coefficients of hydration expansion as the 
most important physical parameters are determined from liter- 
ature data [19]. Strains induced by the hydration change are 
validated using a three-sided constrained membrane in a frame. 
Based on the obtained material properties simulations have been 
run for a four-sided clamped membrane exposed to hydration 
and temperature cycles. Stress—strain results are also compared 
for a clamped membrane with a center hole and results are shown 
for 10 hydration and temperature cycles. In the last section of the 
paper, the time (strain-rate) dependent visco-plastic behavior of 
the membrane is introduced and discussed. 


2. Material properties as a function of T and RH for 
ionomer membrane (Nafion® 111) 


Mechanical properties of Nafion® (in this study commercial 
Nafion® NR111 has been used) are strong functions of relative 
humidity (water uptake) and temperature. Thus, in order to accu- 
rately characterize Nafion®’s properties it is critical to study its 
mechanical behavior for different environments. Uniaxial tensile 
testing was used to investigate the material mechanical prop- 
erties as a function of T and RH. Tests were carried out in an 
environmental chamber (gas chamber) with independent RH and 
T control, built on an MTS Tytron machine. A water bath tray 
was also designed to enable tensile testing with the membrane 
sample fully immersed in water. The specimen was cut into rect- 
angular shape with about 50 mm length and 6.5 mm width. The 
width of membrane sample was measured at three different loca- 
tions across the length of the specimen before the test in room 
condition. All tests were performed at a constant strain rate of 
approximately 0.014s~', but slight changes in strain rate did 
not show significant influence on mechanical properties of the 
material. The original cross-sectional area and gauge length of 
the specimens were used in calculating the stress—strain curves 
from the displacement—load data. 

Fig. 1 shows the engineering stress-strain curves of as- 
received Nafion® NR111 membrane at 25°C with different 
RH from dry to wet (fully immersed in liquid water). Typical 
stress—strain curves of the membrane in water vapor environ- 
ment have a “bilinear” shape consisting of a linear elastic region 
and a linear plastic region, which are interconnected by the 
yielding region, while typical stress—strain curves of membranes 
immersed in liquid water appear to have a power law shape 
with no distinct yield region. The elastic modulus was obtained 
by calculating the initial slope of the stress-strain curve. The 
yield point was defined to be the intersection of extended linear 
parts of the elastic and plastic region. This behavior transition is 
believed to be the result of internal structural change due to water 
uptake. With the increase in water content, there is a decrease 
of both elastic modulus and yield stress. The slope of the plastic 
region (secondary slope) also seems to decrease with increas- 
ing RH. Water, absorbed by polymer, has been recognized as a 
plasticizer, which will soften the material and reduce the load 
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Fig. 1. Stress—strain curves of Nafion® NR111 membrane in 25 °C with varying 
RH. 


carrying capability [48]. Stress—strain curves of the membrane 
in 50% RH with different temperatures from 25 °C to 80°C are 
presented in Fig. 2. It can be seen that the elastic modulus and 
yield stress are decreasing with increasing temperature. Tem- 
perature seems to have a more significant effect on mechanical 
properties of Nafion® membrane than hydration, especially on 
the yield stress. Slight change in the secondary slope is also 
observed when T varies. However, comparing to elastic modu- 
lus and yield stress, the influence of T and RH on the secondary 
slope is less significant. 

Poisson’s ratio was obtained experimentally by using a video 
extensometer and tracing the forward and lateral displacement 
of pre-drawn dots on the membrane during the uniaxial test. It 
should be mentioned that Poisson’s ratio could also be a function 
of hydration and temperature as well. While this issue will be 
addressed in another paper, initial results of our experiments 
have shown that this dependence is not significant and therefore 
it has been assumed to be a constant equal to 0.4 in the present 
work. 
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Fig. 2. Stress-strain curves of Nafion® NR111 membrane in 50% relative 
humidity and various temperatures. 
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3. Incorporation of the material model into FEM 


In PEM fuel cells, the membrane (as part of the membrane 
electrode assembly) is clamped on four edges. This kind of in- 
plane geometry will be addressed in this paper and out-of-plane 
pressure will be ignored. Clamped boundary conditions can be 
the clamped seal regions or the local clamps applied by ribs of 
bipolar plates (flow fields) on the membrane. Therefore, ignor- 
ing the out-of-plane forces does not affect the generality of the 
model. Itis important to point out that there is no external applied 
mechanical force and that changes in hydration and temperature 
are the only driving forces in this simulation. In order to calculate 
accurate numbers for hydration and temperature induced strains 
it is needed to define the coefficients of hydration and temper- 
ature expansion referred to as a and £ in this paper. Having 
proper numbers for alpha and beta, correct values for hygro- 
thermal stresses, which are the major mechanical forces in the 
model, will be calculated by a finite element code according to 
the following expressions in matrix notation [49]: 


{o} = [ST de} — {a} AT — {B}Aa) a) 


where {o} and {£} are mechanical stress and strain tensors due 
to change in hydration and temperature and [S] is the mate- 
rial compliance matrix. œ and £ are the coefficients of thermal 
and hydration expansion and AT and Ad are changes in tem- 
perature and hydration, respectively. Note that if material is 
unrestrained during the hygro-thermal exposure, there are no 
stresses generated and the strains are given by 


{£} = {eng} = {a}AT + {B}Ar (2) 


If the material is completely restrained during hygro-thermal 
exposure, however, the total strain must be zero. Thus, 


{e} = 0 = [S]{o} + {a} AT + {BAA (3) 
and the resulting hygro-thermal stresses are given by: 


{B}AA) (4) 


{o} = [S]"'(-{a}AT 


Establishing the proper value for the coefficient of hydra- 
tion expansion (£) is a key component in the finite element 
model. Nafion® expands/shrinks to a great extent when exposed 
to hydration cycles. Shrink tension and membrane expansion 
induce stress cycles in the membrane during fuel cell operation. 
Changes in dimension with hydration are orders of magnitude 
bigger than we would expect from changes in temperature. In 
the fuel cell this change in hydration (relative humidity) is equiv- 
alent to mechanical load and is believed to be a major driving 
force for mechanical failure. Numbers for coefficients of hydra- 
tion expansion for Nafion® have been extracted from the graphs 
presented by Morris and Sun [19] as shown in Fig. 3. They have 
obtained two curves for change in diameter and thickness of 
membranes. Assuming these two curves are linear, slopes of the 
two lines will be the fs along the diameter and thickness in 
the model, Bx and By. These two values are 0.013 and 0.008, 
respectively, and have been incorporated into the model — along 
with coefficients of thermal expansion a — by the ABAQUS 
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Fig. 3. Coefficients of hydration expansion along the diameter and thickness of 
Nafion®. X-axis is mass ratio of water to Nafion® fm, while Y axis is dimensional 
change. Water content of membrane (A) is equal to 61.0 x fm [19]. 


user subroutine UEXPAN [50]. For all the in-plane simulations 
presented in this paper just the first value is used. 

Material properties (uniaxial stress-strain curves) for the 
Nafion® as shown in the previous section are obtained as a 
function of temperature and relative humidity. Hence, water con- 
tent, as an independent field variable in finite element modeling, 
must be related to relative humidity for different temperatures. 
Springer et al. found a relationship between water content and 
relative humidity of N afion® [51]. This equation (Eq. (5)) relates 
relative humidity (water vapor activity) to water uptake from 
water vapor at room temperature (30°C). 

There is controversy in the literature on the relationship 
between Nafion®’s water content and relative humidity with 
change in temperature. Some researchers found that for a cer- 
tain relative humidity, water uptake decreases when temperature 
increases [18,52—54] while others reported the opposite trend 
in which water uptake goes up as temperature rises [19,55,56]. 
Hinatsu et al. reported the relationship for different temperatures 
pointing out that the shape is similar for all the temperatures [18]. 
Eq. (6) is the equation they derived for T= 80°C which shows 
that water content at 80 °C is less than that reported for the lower 
temperatures. They have related this temperature effect to diffi- 
culty in condensing water on the membrane surfaces at higher 
temperatures. Eqs. (5) and (6) will be used in this paper to relate 
à (water content) to relative humidity (water activity). 


@T = 25°C: à = 0.043 + 17.81a — 39.85a7 + 36.0a> (5) 
@T = 80°C: A = 0.300 + 10.8a — 16.0a? + 14.14? (6) 


Values for other temperatures can be interpolated (extrapo- 
lated) from these two equations. The value of the coefficient of 
thermal expansion (œ) is chosen to be 1.47 x 1074 [57]. 

These equations have been obtained for membranes in contact 
with water vapor, which is the primary water uptake mechanism 
for fuel cell humidification. Water uptake from the liquid phase 
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with unit water activity (100% relative humidity) is higher than 
water uptake from the saturated vapor phase, a phenomenon 
known as Schroeder’s paradox [58]. This phenomenon is the 
main reason why water content is used as a field variable instead 
of relative humidity in the numerical simulations. This pre- 
vents ambiguities when it comes to fully saturated environments 
(RH = 100%, vapor or liquid phase) and numerical convergence 
problems with associated computations are unlikely to happen 
for that reason. For instance, for unit water activity at room tem- 
perature, the water content in the membrane can be any number 
from 14 to 22; this imposes different swelling conditions and 
therefore different stress states in the membrane, which is hard 
to represent in the numerical code. It has been reported in liter- 
ature that unlike the case of water uptake from the vapor phase, 
water uptake in the liquid phase goes up with increasing tem- 
perature. Membranes placed in water at 100°C for 1 h absorbed 
almost twice the amount of water as membranes placed in water 
at 25 °C for 24 h [24]. This issue must be addressed in water elec- 
trolyzers where the membrane is in contact with liquid water, but 
will not be considered here. Physically, a membrane sandwiched 
between two bipolar plates has a lower value of à compared to 
a free swollen one that is not constrained. Pressure from bipolar 
plates in a fuel cell configuration can push some water out of the 
membrane in a similar way as a squeezed sponge [44,47]. This 
effect has not been accounted for in the presented simulations. 

In order to model the strains/stresses developed during the 
RH and temperature cycling of a clamped membrane, hydration- 
dependent changes of material properties during RH (and 7) 
variation have been incorporated into the commercial finite ele- 
ment (FE) code ABAQUS (HKS). The measured T and RH (or 
à) dependent stress-strain responses discussed in the previous 
section of this paper were incorporated into the simulation model 
through the ABAQUS built-in elasto-plastic material model. The 
4-dependent material properties were implemented through look 
up tables, such as: 


*ELASTIC *PLASTIC 
E1, vı, Ti, M41 O1, Ep1; Ti, À1 
Fo, v2, T2, d2 02, Ep2, T2, A2 


where E and v=0.4 are the elastic modulus and Poisson’s ratio, 
respectively; o and £p are true stress and strain extracted from 
the plastic region of the stress-strain curve; T and A are the cor- 
responding temperature and water content (RH). The following 
equations are used to convert engineering stress-strain to the 
true stress—strain values needed in FEA. 


Etue = ln(1 + Eengr) (7) 


Otrue = Sengr oe? In(1+ Seng) (8) 


As the applied RH changes, ABAQUS updates the materials 
constitutive behavior according to the table. There is no exter- 
nal mechanical force applied in the simulation of the clamped 
membrane and the only driving force is the change in hydration 
as an independent field variable. It can be easily concluded from 
the values of œ and 6 and Eqs. (1) and (2) that hydration has 


a far bigger effect than temperature on shrinkage/expansion of 
membranes in the fuel cell environment. For this reason, change 
in hydration will be the primary driving force in all the present 
simulations. At the same time, change in temperature will also 
be applied in order to pick the appropriate material curves dur- 
ing the simulation. It is also known that Nafion®, as a polymer, 
demonstrates time-dependent behavior such as creep and stress 
relaxation. That behavior was not incorporated in this reported 
model, but it can be included in the model and will be discussed 
in a future work by adding a viscous module to this elasto-plastic 
model. 


4. Validation of strains and hydration coefficients 


In the experiment conducted to verify the FEM results, a rect- 
angular fully hydrated specimen was clamped in a rigid frame. A 
26mm x 23 mm specimen of NR111 was used. Three edges of 
the specimen were constrained and black dots have been made on 
one edge of the specimen to be traced by the video-extensometer. 
After waiting for the dots to dry, the specimen was hot pressed at 
100°C to ensure the dots adhered to the membrane. The mem- 
brane was then boiled in DI water for about an hour. It was 
then placed in the designed frame such that three of the mem- 
brane edges were constrained and the dotted edge remained free. 
During installation of the specimen on the frame, two pieces of 
wet Kimwipe™ were placed on both sides of the membrane to 
reduce water loss during mounting. After tightening the bolts 
of the frame the Kimwipes™ were removed and extra water on 
the surface of membrane and also on the frame was blotted with 
Kimwipe™. The front surface of the clamped membrane was 
viewed by a video-extensometer to trace the displacements of 
the dots on the free edge with time while the membrane dried 
out from fully hydrated to room condition. It can be seen that the 
free edge of the membrane shrinks and deforms into a meniscus 
shape as it loses water. Fig. 4 shows the constrained membrane 
in the frame at the end of the experiment. 

A rectangular specimen of a commercial ionomer membrane 
with three constrained edges and one free edge (on top in the 
figure) was modeled with ABAQUS. Material curves for 23 °C 
and hydrations of 25%, 50%, 75% and 100% (underwater tests) 


Fig. 4. Clamped membrane in the frame at the end of test. 
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Fig. 5. Comparison of model prediction and experimental measurements. 


were incorporated into the model. The applied ‘load’ was the 
change in hydration, from 100% to 60% in this case, which was 
chosen based on the initial and final conditions of the membrane 
in the validation experiment that was conducted. Fig. 5 shows 
the deformed meniscus from the experiment along with the pre- 
dictions of the FEM models. As seen in the picture, by using 
the proper values for hydration expansion coefficient and mate- 
rial properties correct deformation states of the material can be 
predicted by using the finite element method as discussed above. 


5. Simulation results for a constrained membrane 
subjected to change in hydration and temperature 


It is the purpose of the simulations presented in this section 
to establish the shrink/expansion stresses and their distributions 
in a framed membrane subjected to changes in hydration and 
temperature. As mentioned before, in the fuel cell configuration 
mechanical failure occurs without external mechanical loads 
applied in the plane of membrane. The only applied mechan- 
ical load is the out-of-plane clamping pressure from the bipolar 
plates. Hence, the only forces developed in in-plane directions 
are the tensile/compressive forces due to membrane shrink- 
age/swelling. However, there is always a contribution from 
the out-of-plane clamping pressure as a form of compressive 
stress-strain in the plane of membrane. Out-of-plane pressure 
tries to squeeze the membrane and therefore pushes the mem- 
brane towards the clamped walls. Reaction forces from the walls 
result in some compressive stress in the plane of the mem- 
brane which will be ignored in the simulations. Plane stress 
simulations are performed. Therefore, there is no out-of-plane 
contribution in the model and resulting stresses are solely due to 
in-plane swelling/shrinking behavior. Fig. 8 shows the geom- 
etry (30cm x 30cm square) and boundary conditions of the 
finite element model. This model consisted of 400 linear plane 
stress elements (441 nodes) with clamped boundary conditions 
at the edges. Clamped boundary conditions are not necessar- 
ily the seal locations in the fuel cell assembly. The clamped 
edges could be a result of any applied force out of plane of the 
membrane that locally constrains parts of the membrane. These 


locally constrained geometries could be the results of out-of- 
plane forces from ribs of bi-polar plates over the gas diffusion 
layer and, therefore, the membrane. Applied loads are changes 
in water content (A) and temperature defined as independent 
field variables. Elasto-plastic material models are incorporated 
as explained in Section 3. This elasto-plastic model does not 
consider the time (rate) response of the polymer membrane; a 
visco-plastic model is needed to capture time dependency of the 
material response and will be added subsequent to this work. 
Hence, the model physically corresponds to a relatively fast dry- 
ing process close to the uniaxial test strain rates. However, the 
physical uptake or egress of water in Nafion is quite rapid. 

Fig. 6 shows the stress results as temperature and hydration 
changes are applied to the membrane. The authors have reported 
previously that crazing is one of the mechanisms contributing to 
failure of the membrane electrode assemblies [28]. Since crazing 
is initiated by hydrostatic stress [59] it is essential to study the 
change in this variable as well as in the deviatoric term as a 
function of temperature and hydration. It has to be mentioned 
that since all the analyses in this paper are plane stress and hence 
033 is assumed to be zero, the results reported in this paper 
for stresses, especially hydrostatic stress (oj;) are conservative 
results and do not capture all of the clamping pressure effects of 
the bi-polar plates. In the case shown in Fig. 6, the water content 
is changing from | to 9.2 (water content at fully saturated vapor 
phase for 80°C temperature calculated from Eq. (6)) while the 
temperature changes from 25 °C to 80°C. Stress-free condition 
is the dry state where temperature equals 25 °C and water content 
value is 1. This is analogous to the in-plane condition in which 
the fuel cell starts up from a relatively dry condition at room 
temperature to a saturated (wet) mode at operating temperature. 
It is clear that swelling stress is building up in the membrane 
with increase in water hydration as temperature also goes up. 
At some point during hydration, the stress relaxes to a lower 
value. This could be as a result of the lower stiffness material 
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Fig. 6. Von-Mises and hydrostatic stress developed in the membrane as a result 
of change in hydration (A: 1 > 9.2) and temperature (T: 25°C > 80°C). 
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Fig. 7. Hydrostatic stress for different cases of temperature and hydration swing. 


properties at higher temperature and moisture as mentioned in 
Section 2. Basically, smaller stiffness values compensate some 
of the stress induced by hygro-thermal expansion. Formation of 
permanent (plastic) deformations due to the increase in hygro- 
thermal strains also contributes to this relaxation of stress with 
increases in A and T. 

Obviously, since the membrane is constrained on four sides 
and there is no non-uniformity in the plane of the membrane, 
there will not be any resulting in-plane total strain in the model. 
However, because of the positive values for hygro-thermal 
strains as a result of expansion tendency of the membrane, 
there will be negative elastic and plastic strains in the model 
to maintain zero total strain. Compressive plastic strain during 
hygro-thermal loadings has been reported in literature [30]. 


{Etotal} = 0 = {Eelastic} + {€plastic} + {€hygrothermal } (9) 


Since the material model at this point does not account for 
compression, these negative strains will not be further discussed. 
Fig. 7 compares the simulation results when just temperature 
or hydration or both are changing. This illustrates the effects 
of hydration and temperature swings independently in the sim- 


Outlet 


30 cm 


ulation. In the temperature-only case A is fixed at 1 and the 
temperature is changing from 25°C to 80°C, while in the 
hydration-only case, the temperature is fixed at 25°C and the 
hydration swing is from 1 to 9.2. It is clear that cycling à has 
a far bigger effect than temperature cycles. These results reveal 
the counter-intuitive nature of the problem when A, T and mate- 
rial properties change at the same time. From Fig. 7 it can be 
seen that stress, with just increase in T or A, is monotonically 
increasing, but it does not follow the same trend indicated when 
both T and å increase simultaneously. This observation points 
out that simple superposition of each individual case with just 
increasing T or 4 can not be used to get the correct results when 
both of them increase or decrease at the same time, and it shows 
the significance of the finite element analysis as a method to 
obtain the correct behavior. 

For all these cases, the simulation results are shown for a 
node in the middle of the membrane. The simulation results are 
uniform in the plane of the membrane since there is no non- 
uniformity in the applied conditions. Generally, in a membrane 
operating in a fuel cell, point-to-point variations in temperature 
and hydration of the membrane are inevitable. The membrane 
tends to be more greatly humidified near the outlet of the flow 
channels as a result of water production when fuel cell reactions 
take place. 

This non-uniformity dictates different point-wise stress— 
strain behavior in the membrane during operation. In order to see 
this effect, a hydration profile is applied across the membrane 
as shown in Fig. 8. In this case same temperature and hydration 
change as discussed above was applied, but at the end of these 
changes hydration was increased near the end of the flow chan- 
nel (top edge in Fig. 8) from 9.2 to 20 to simulate the possibility 
of the liquid phase water increasing at the outlet. Stress-free 
condition is again T7=25 °C and à = 1. Figs. 9 and 10 show the 
stress and strain contours of the membrane as a result of this 
hydration profile. The maximum values for total strain are near 
the more hydrated edge of the membrane. The strain developed 
in the plane of the membrane in the constrained configuration 
could raise wrinkling issues in the membrane near the edges 
during fuel cell operation. Since material cannot expand in the 
plane of the membrane, it is possible that it would wrinkle out 
of its plane due to hydration expansion. Unlike the uniform case 
where no shear effect was observed, non-uniform hydration also 
develops shear stress (as well as shear strains) with maximum 
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Fig. 8. Generated mesh and schematic of the applied hydration profile as a result of water generation towards the outlet of flow channels. 
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Fig. 9. Contour plot of hydrostatic stress (a) and maximum in-plane principal strain (b) when a linear hydration profile is applied across the membrane. Maximum 
value for hydrostatic stress is 4.613 MPa near the corners and maximum principal strain is 0.19 around the middle of top edge. 


values near the left and right edges (Fig. 10). This observed shear 
would not be expected for a uniform symmetric geometry. Max- 
imum values of hydrostatic stress happen near the top corners 
and can be considered as potential spots for failure initiation 
(Fig. 10a). These results clearly show the significance of the 
non-uniformities across the membrane, which have to be con- 
sidered for further understanding of the membrane’s mechanical 
behavior and potential failure modes. 

Creation of a hole in the membrane can possibly be consid- 
ered as the extreme form of non-uniformity in this constraint 
configuration. A hole in the membrane can represent pinhole 
creation or even material degradation during fuel cell operation. 
As part of this modeling work a membrane with a circular hole in 
the center has been modeled and stress-strain results are studied 
as a function of change in environmental variables. Geometry is 
again 30 cm x 30 cm square but this time with a 1.5 cm diameter 
hole in the center and boundary conditions are again clamped 
around the edges. Generated mesh has 3608 nodes and 3504 
plane stress elements with finer mesh near the hole to capture 
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stress-strain concentrations. Fig. 11 shows the geometry, mesh, 
and boundary conditions applied to the membrane. 

Two different cases will be studied for the center-hole con- 
figuration. Firstly, we will consider the uniform change from the 
stress-free state at T=25 °C and A= 1 to T=80°C and à =9.2 as 
it was done before. In the second case, after the first step (mod- 
eled in the first case), environmental conditions change in the 
second step from the previous hot, hydrated state (T= 80 °C and 
à =9.2) to the initial cold and dry state which was T= 25 °C and 
à = 1. These conditions simulate the start-up and shut-down of an 
operating fuel cell. In the second step the material shrinks while 
in the first step the material tends to expand. These two steps 
are repeated five times (equivalent to 10 total start-up and shut- 
down cycles). This process estimates the accumulative plastic 
strain developed in each step which can be studied. 

Fig. 12 shows the distribution of stress and strain in the mem- 
brane with the center hole. From Fig. 12, it can be seen that the 
stress concentrations are 1.38 for Mises stress and 0.7 for hydro- 
static stress component. It should be noted that in the elastic case 
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Fig. 10. Contour of shear stress (a) and shear strain (b) developed as a result of non-uniformity in moisture distribution across the plane of membrane. Maximum 
value for shear stress is 0.3 MPa and for shear strain is 0.02 both near the middle of side edges. 
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30 cm 


Fig. 11. Geometry and mesh of the constrained membrane with a center hole. 


the Mises stress concentration would be 2 and there would be no 
variation in hydrostatic stress. This can be easily obtained from 
the elasticity solutions of a biaxially loaded plate. Fig. 12b shows 
the distribution of maximum principal strains along the center- 
line of the membrane. As mentioned before, the farfield values 
for plastic and elastic strains are negative because of the con- 
stant positive value for hygro-thermal strain and the constrained 
geometry (Eq. (9)). 

Stress and strain (maximum principal) values at a single node 
at the edge of the hole (as indicated in Fig. 11) for 10 loading 
and unloading cycles are shown in Fig. 13. During the unloading 
cycle (shut down), while the membrane expansion in the load- 
ing cycle is shrinking back to the initial dry state, von-Mises 
stress reaches a minimum value and then at some point increases 
to a maximum value at the end of the cycle. This maximum 
value is relatively high and can exceed the ultimate stress values 
for Nafion® (the maximum von-Mises stress is about 22 MPa 
in this simulation). Ultimate stress and strain-to-failure values 
will be reported and discussed further in authors’ other publi- 
cations. This unsymmetrical behavior in loading and unloading 
is because of the nature of unrecoverable strains developed in 
the loading cycle. Similar asymmetry is observed in hydrostatic 
stress. 
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The most interesting observation is in the plastic strain results. 
After the first three load/unload cycles plastic strain at the edge 
of the hole reaches a constant value and does not fluctuate with 
changes in hydration and temperature. It suggests that there is a 
saturation value for plastic strain which occurs after a few cycles. 
This can be associated with the uniform permanent change in 
the shape of the circular hole as a result of plastic deformation. 
It can also be seen that other strain and stress results will show 
similar fluctuations after the first six cycles (three load/unload) 
cycles because of the constant plastic strains. 


6. Visco-plastic model for membrane 


Nafion® as a polymer shows strong time dependency. Fig. 14 
shows the uniaxial mechanical response of Nafion® at room 
condition (T=25 °C and RH=50%) for different initial strain 
rates. In each of these curves the specimen is mounted and equi- 
librated in the environmental chamber at 25 °C and 50% RH for 
1 h before the test is performed. The specimens were pulled with 
different initial strain rates up to an engineering strain of 0.5 and 
then the displacement was fixed for 1 h in order to observe the 
stress relaxation. It is clear from this figure that significant over- 
stress is observed for different strain rates and the stress value 
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Fig. 12. Stress (a) and strain (b) distribution in the membrane with a hole in the center. X-axis shows the distance from the hole in the centerline of the membrane. 


Strain values are maximum principal strains. 
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Fig. 13. Stress and strain results at the edge of the hole after ten total load and unload cycles. Plastic strain will be saturated to a constant value after about six cycles. 


Strain values are maximum principal strains. 


relaxes to an equilibrium state with time. Therefore, a classi- 
cal elasto-plastic model, i.e., classical isotropic plasticity, is not 
sufficient to capture the accurate stress values when the material 
exhibits time and rate effects. Stress relaxation for one of the 
four curves is depicted in Fig. 15. 
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Fig. 14. Stress—strain curves of Nafion® at room conditions for different strain 
rates. 


Previous work reported on Nafion®’s viscous behavior has 
been reported in the literature in which a linear visco-elastic 
model, which is accurate in the linear elastic region when applied 
strains are relatively small, has been used [29]. By assuming 
visco-elastic behavior for Nafion®, modeling its behavior cannot 
be accurately done when the material undergoes large strains and 
is therefore viscoplastic. Because of the large swelling charac- 
teristic of membranes in fuel cell/electrolyzer applications, large 
strains and, therefore, considerable unrecoverable deformations 
are inevitable. 
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Fig. 15. Stress relaxation of Nafion® at room conditions. Initial applied true 
strain is 0.37 with strain rate of 2.34e-2. 
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The authors realize that amore complete viscoplastic model is 
required to predict stress-strain values for the investigated poly- 
mer membranes. In a future paper a two-layer viscoplasticity 
model will be proposed based on the observed time-dependent 
behavior of the material and availability of the model in the com- 
mercial FEA package ABAQUS (HKS) be used in conjunction 
with the elastic—plastic module developed in this work. Simula- 
tion results based on this model show good agreement with the 
results shown in Figs. 14 and 15. 


7. Conclusion 


By incorporating the non-linear material curves for Nafion® 
NR111 obtained from controlled environment uniaxial exper- 
iments in the laboratory and coefficients of hygro-thermal 
expansion from the literature, an elasto-plastic model was 
developed and used to estimate the mechanical behavior of 
Nafion® membrane using the commercial finite element package 
ABAQUS. A non-linear visco-plastic model is also needed for 
precise time-dependent modeling of ionomer membranes. How- 
ever, elasto-plastic simulations still provide valuable insights for 
understanding the RH-induced straining of the membrane in the 
constrained fuel cell configurations. 

Model simulations show that hydration and drying cycles due 
to start-up/shut-down of fuel cells can cause significant stresses 
in the clamped membrane. It is essential to study these induced 
mechanical stresses to understand the failure of the clamped 
membrane, where there is no significant external mechanical 
loads. Hydration is shown to have a bigger effect than tempera- 
ture in developing mechanical stresses in the membrane. These 
stresses will be more critical when non-uniformity as a form of 
hydration profile or a physical pinhole exists across the mem- 
brane. As a result of a hydration profile (due to water production 
in the flow channels), non-uniform stress and strain results were 
observed in the plane of the membrane with maximum values 
near the edges. Because of this non-uniformity, shear is seen in 
the plane of the membrane even though for a symmetric uniform 
geometry shear effects are usually not expected. 

Values for stress and strain concentrations have been obtained 
near the circular hole formed in the membrane when subjected to 
changes in hydration and temperature. Because of the permanent 
deformations (plastic strain), concentration values for deviatoric 
and hydrostatic stresses are lower than the elastic case. How- 
ever, in the case of a discontinuity (pinhole) in the plane of the 
membrane, stress values can exceed the ultimate stress of the 
material and it is likely that a crack initiates at the edge of the 
pinhole. By repeating the hydration/dehydration cycles, plastic 
strain will saturate to a constant value and does not fluctuate 
after about three load/unloading (six total) cycles in the present 
case. 

Stress—strain results of the constraint configuration due to 
hydration/dehydration discussed in this paper provide insights 
to better understand the behavior of ionomer membranes that are 
clamped and hydrated in fuel cells. These results show the impor- 
tance of hydration swelling and shrinkage when a new polymer 
membrane is designed. It is crucial to design and control the 
expansion properties of a membrane in addition to its electro- 


chemical properties in the manufacturing process. By reducing 
the hygro-thermal dimensional change of a membrane, mechan- 
ical stress/strains can be reduced which can extend the life of 
the membrane. Studying in-plane mechanical response of mem- 
brane as well as through-the-thickness stress—strain behavior 
and failure analysis of membranes in the membrane electrode 
assembly due to mechanical cycling is essential to investigate 
mechanical degradation of membranes in fuel cells. Detecting 
the driving forces for mechanical degradation and also chemi- 
cal degradation of the membrane is necessary to determine the 
relevant failure modes of such membranes. Mechanistic based 
failure metrics based on experimental and numerical observa- 
tions and the present modeling approach will be used in the 
authors’ future work for damage accumulation analysis in dura- 
bility modeling and life prediction of fuel cells. 
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